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1. Introduction 



First data reported by the STAR Collaboration at RHIC P has a signif- 
^ I leant meaning that the observed large magnitude of elliptic flow for charged 

Ch ' hadrons is consistent with hydrodynamic predictions i 2_ . This suggests that 

large pressure possibly in the partonic phase is built at the early stage 
(r ~ 0.6 fm/c) in Au+Au collisions at ^/s^ = 130 and 200 GeV. This 
situation at RHIC is in contrast to that at lower energies such as AGS or 
I SPS where hydrodynamics always overpredicts the data [H]. Moreover, this 

also suggests that the effect of the viscosity in the QGP phase is remarkably 
small and that the QGP is almost a perfect fluid [4 . Hadronic transport 
models are very good to describe experimental data at lower energies, while 
they fail to reproduce such large values of elliptic flow parameter at RHIC 
(see, e.g., Ref. j^]). So the importance of hydrodynamics is rising in heavy 
ion physics. After the first STAR data were published jlj, other groups at 
RHIC have also obtained the data concerning with flow phenomena [H]. 

Contrary to the success of the hydrodynamics to describe the elliptic 
flow, many dynamical models including hydrodynamics cannot reproduce 
the HBT radh [31311]. It is known as the "HBT puzzle". 
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To understand these experimental data, hydrodynamic analyses are per- 
formed extensively IQ, ■ this short review, we highlight several results 
mainly on elliptic flow and on the HBT radii from hydrodynamic calcula- 
tions. 

2. Basics of Ideal Hydrodynamics 

Hydrodynamic equations represent the energy-momentum conservations 
d^T^^'^ = 0. In the ideal hydrodynamics, the energy-momentum tensor 
becomes T^'^ = (e + P)u^vy — Pg^" , where e is the energy density, P is 
the pressure, and u'^ is the four fluid velocity. When there are conserved 
quantities such as the baryon number or the number of chemically frozen 
hadrons, one needs to solve the continuity equations 5^nf = together 
with the hydrodynamic equations. In order to close the system of partial 
differential equations, the equation of state (EoS) P(e, rij) is needed. The 
naive applicability conditions of ideal hydrodynamics are that the mean free 
path among the particles is much smaller than the typical size of the system 
and that the system keeps local thermal equilibrium during expansion. From 
these conditions, one cannot use hydrodynamics for initial collisions, final 
free streaming and high pT 2 GeV/c) particles. So one needs an interface 
between the pre-thermalisation stage and the hydrodynamic stage at the 
initial time. Moreover, the system eventually breaks up and cannot keep 
thermalisation at the later stage. This means a prescription of freezeout 
is needed in the hydrodynamic model in evaluating the particle spectra. 
Therefore one needs another interface between the hydrodynamic stage and 
the free streaming stage. In what follows, we discuss particularly equations 
of state, initial conditions and freezeout prescriptions used in the literature. 

2.1. Equation of State 

The main ingredient of the hydrodynamic model is the equation of state 
(EoS) for thermalised matter produced in heavy ion collisions. Ideally, one 
uses the EoS taken from the first principle calculations of QCD, namely, lat- 
tice QCD simulations 1^. More practically, one can use the resonance gas 
model for the hadron phase and the massless free parton gas for the QGP 
phase. By matching these two models at the critical temperature, one ob- 
tains the first order phase transition model with a latent heat ~1 GeV/fm^. 
In the mixed phase at = 0, the sound velocity Cg = dP jde is vanishing. 
Recent lattice QCD simulations tell us that the phase transition seems to 
be crossover in vanishing baryonic chemical potential. Discontinuity of the 
thermodynamic variables does not exist in the crossover phase transition. 
It should be emphasized, however, that the energy density and the entropy 
density increase more rapidly than the pressure in the vicinity of the phase 
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transition region AT ~ O.lTc. This also leads to very small sound velocity 
near the phase transition region. Therefore it is very hard in general to 
find flow observables which distinguish the crossover phase transition with 
a rapid change of the thermodynamic variables from the first order phase 
transition. 

2.2. Initial Condition 

Once initial conditions are assigned, one can numerically simulate the 
space-time evolution of thermalised matter which is governed by hydrody- 
namic equations. Usually, initial conditions are parametrised based on some 
physical assumptions. Transverse profile of the energy /entropy density is 
assumed to be proportional to the number density of participants Ppart, 
the number density of binary collisions p^ow or linear combination of them. 
Initial transverse flow is usually taken to be vanishing. 

On the other hand, one can introduce model calculations to obtain the 
initial condition of hydrodynamic simulations. Event generators can be used 
to obtain the energy density distribution at the initial time. Recently, the 
SPheRIO group employs an event generator NeXus and takes an initial con- 
dition from this model in the event-by-event basis |13[ I14j . The resultant 
energy density distribution in the transverse plane has highly bumpy struc- 
tures jl31 115j . Smooth initial conditions used in the conventional hydrody- 
namic simulations are no longer expected in one event. Another important 
example which is relevant at very high collision energies is an initial con- 
dition taken from the Colour Glass Condensate picture. See Ref. [IB] for 
recent calculations. 

2.3. Freezeout 

Conventional prescription to obtain the invariant momentum spectra 
from the hydrodynamic simulations is to employ the Cooper-Frye formula 
jl7j . The physical picture described by the Cooper-Frye formula is some- 
times called "sudden freezeout" since the mean free path is suddenly changed 
from zero to infinity through a thin freezeout hyper surf ace. Instead of using 
this, one can use a hadronic cascade model to describe the space-time evo- 
lution of hadrons jl8l I19j . The mean free path among hadrons is finite and 
depends on hadronic species. Hence, one can describe a continuous freezeout 
picture through hadronic transport models. Note that continuous particle 
emission can be considered within the hydrodynamics ^20i. Adoption of 
hadronic transport models after hydrodynamic evolution of the QGP liq- 
uid could refine the dynamical modeling of relativistic heavy ion collisions. 
However, it is not so easy to connect them in a systematic and proper way 
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Fig. 1. (Left) Transverse momentum dependence of V2 for pions and (anti)protons 
[TT] are compared with the STAR data in minimum bias Au+Au collisions at 
^/SNN = 130 GeV. EoS Q stands for the first order phase transition model, while 
EoS H stands for the resonance gas model. (Right) V2 for charged hadrons as a 
function of pseudorapidity rj in Au+Au collisions at y/sWN = 130 GeV [Tl 1^. 
PCE means the EoS of partial chemical equilibrium, whereas CE means the EoS 
of chemical equilibrium. Figure is taken from Ref. j24| . 

3. Hydrodynamic Results for V2 

Assuming the Bjorken flow [22. for the longitudinal direction, one can 
solve the hydrodynamic equations only in the transverse plane at midra- 
pidity. Systematic studies based on this (2+l)-dimensional hydrodynamic 
model are performed in Ref. |2j. For the EoS, complete chemical equilibrium 
is assumed for both the QGP phase and the hadron phase. pT dependences 
of V2 for pions and protons from this model Tl, are compared with the STAR 
data [2'6\ in Fig. ^ (left). By employing the EoS with phase transition, the 
hydrodynamic model correctly reproduces V2{pt) and its mass-splitting be- 
havior below pt = I GeV/c. On the other hand, V2{pt) for (anti)protons 
from the resonance gas model does not agree with the data. Although the 
reason for the difference of the result between these two EoS models is not so 
clear, the experimental data favors the QGP EoS. Due to the assumption 
of chemical equilibrium in the hadron phase, this model does not repro- 
duce particle ratio and spectra simultaneously. It is of importance to study 
whether the agreement with the experimental data still holds even when the 
assumption of chemical equilibrium in the hadron phase is abandoned . 

One needs a full 3D hydrodynamic simulation to obtain the rapidity 
dependence of V2. First analysis of V2{ti) at RHIC based on the full 3D 
hydrodynamic model is performed in Ref. i211 [25j . Figure Q (right) shows 
V2{ri) for charged hadrons in Au+Au collisions at -^Jsnn = 130 GeV [TIOS]. 
Here the initial condition for the energy density is so chosen as to reproduce 
the pseudorapidity distribution of charged hadrons. Hydrodynamic results 
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Fig. 2. (Left) Excitation function of V2/S. Figure taken from Ref. |271. (Right) 
Excitation function from the hydro+cascade model. Figure taken from Ref. |19| . 



are consistent with the experimental data only near the midrapidity while 
hydrodynamics overpredicts the data in the forward/backward rapidity re- 
gions. Multiplicity is not so large in the forward/backward rapidity regions, 
so equilibration of the system tends to be spoiled. 

Figure 121 (left) shows the excitation function of V2 compiled by the STAR 
Collaboration j27j . Hydrodynamic results presented in this figure are based 
on the same model discussed in Fig. ^ (left). Data points continuously 
increase with the unit rapidity density at the AGS, SPS and RHIC energies. 
However, the hydrodynamic response f2/e is almost flat or slightly decreases 
with the multiplicity. The data points seem to reach the "hydrodynamic 
limit" for the first time at the RHIC energy. The deviation between the 
hydrodynamic results and the data plots below {1/ S)dN / dy = 25 reminds 
us the pseudorapidity dependence of elliptic flow in Fig. ^ (right). The 
deviation might come from a common origin j28j : The small multiplicity 
both in forward rapidity region at the RHIC energy and at midrapidity at 
the lower collision energies could cause the partial thermalisation. In the 
low density limit, V2 is actually proportional to the number density j29j . 
The shape of V2{r]) data in forward rapidity region looks similar to that 
of the pseudorapidity distributions |3D| . Similarly, data plots of excitation 
function increase almost linearly with the particle density. These results 
suggest that thermalisation is not achieved completely in forward rapidity 
region at the RHIC energy and at midrapidity at the SPS energies. 

Figure |21 (right) shows the excitation function from the hydro+cascade 
(RQMD) model JHI with Bjorken longitudinal flow j22j . Contrary to the 
conventional hydrodynamic models, freezeout processes are automatically 



6 



HIRANOISMD04 PRINTED ON FEBRUARY 8, 2008 



described by the cascade model without any adjustable parameters. The 
excitation function in the case of the latent heat ~ 0.8 GeV/fm^ linearly 
increases with the multiplicity, which is consistent with the experimental 
data. The difference from the conventional hydrodynamic results might 
come from the strong viscosity in the hadron phase within the cascade cal- 
culation. It should be noted that V2{pt), its mass dependences and particle 
spectra/ratio at midrapidity are also reproduced by this hybrid model jl9j . 

4. Hydrodynamic Results for HBT Radii 

Although many hydrodynamic calculations are already performed, one 
does not succeed to interpret the HBT puzzle yet. The main reason why 
the hydrodynamic simulations overestimate the -RoutZ-^side is the negative 
value of the correlation between Xout and t which comes from the hydrody- 
namic source function . Here x = x — {x) and the average is taken over 
the source function. Note that positive Xout-t correlation can be obtained 
dynamically in a transport model So detailed comparison of this re- 
sult with hydrodynamic results would be important in understanding the 
space-time evolution of matter in heavy ion collisions. 

As discussed in Sec. 12.31 the conventional hydrodynamics assumes the 
sudden freezeout which must be far from the realistic situation. The HBT 
radii reflect the distribution of the final scattering points. So realistic treat- 
ment of the final decoupling stage is mandatory. ^ To remove this unwanted 
feature in the hydrodynamic model, the continuous particle emission is pro- 
posed in Ref. |141 l2Uj . This prescription naturally gives that larger mo- 
mentum particles comes from the earlier stage and that smaller momentum 
particles comes from the later stage. For details of the results, see Ref. jl4j . 

As shown in the previous section, the hydrodynamics with a hadronic 
cascade model at the late stage reproduces the momentum space of particle 
spectra from SPS to RHIC energies. This indicates the dilute hadronic 
gas should not be described by the ideal hydrodynamics. In Ref. a 
hydro+cascade (UrQMD) approach is employed to calculate the ratio of 
the pion HBT radii as shown in Fig. |31 In the case of Tc = 160 MeV, the 
ratio becomes 2.0 around pair transverse momentum Kt =0.15 GeV/c at 
hadronization which is much larger than the experimental data ~ 1 [7||Hj. 
Although smearing of freezeout hypersurface by using a hadronic cascade 
reduces the ratio to 1.4-1.6, it is not enough to reproduce the experimental 
data. 

^ The effect of final multiple scattering on the HBT radii is recently discussed in 
Refs. |32l 1331 . According to these analyses, the distribution one can obtain from 
the HBT analysis might be the initial effective one rather than the one of final scat- 
tering points. 
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Fig. 3. i?out/.Rside at hadronization (lines) and at freezeout (symbols) from a hy- 
dro+cascade approach. 



The most successful hydrodynamic approach to elliptic flow in relativis- 
tic heavy ion collisions is the hybrid model in which the QGP phase is 
described by the ideal hydrodynamics while the hadron phase is described 
by a hadronic cascade. However, even within this model, the ratio of the 
HBT radii is still larger than the data. The Bjorken scaling solution is 
assumed in the current hydro+cascade models. This means that current 
hybrid models are available only near midrapidity. Therefore it is desired 
to develop a model in which a full 3D hydrodynamic simulations combined 
with a hadronic cascade model. Prom agreement of excitation function be- 
tween the hybrid model and the experimental data at midrapidity, the 3D 
hybrid model is expected to reproduce the pseudorapidity dependence of 
elliptic flow. It should be emphasised again that the hybrid model has its 
own problem on the violation of energy momentum conservations between 
the QGP liquid and the hadron gas. A hybrid simulation which incorpo- 
rates a proper treatment at the boundary between the QGP phase and the 
hadron phase is now an open problem. 
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